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ANGLE  ESTIMATION  IN  THE  PRESENCE 
OF  MAINBEAM  INTERFERENCE 


I.  INTRODUCTION 

The  problem  of  target  tracking  in  the  presence  of  mainbeam  interference  is  of  considerable 
interest  in  radar  systems.  Conventional  tracking  radars  may  experience  serious  errors  in  this  situa¬ 
tion.  This  problem  was  first  addressed  in  the  literature  by  Davis  et  al.  [1],  who  extended  the  theory 
of  adaptive  arrays  to  the  angle  measurement  problem  based  on  maximum  likelihood  theory.  They 
proposed  an  angle-of-arrival  estimator  involving  adaptively  distorted  sum  and  difference  beams  that 
are  analogous  to  those  used  in  conventional  monopulse  antennas.  Good  per  ormance  by  simulation 
was  demonstrated  for  sidelobe  and  mainbeam  interference. 

The  approach  taken  by  Gabriel  [2]  was  to  determine  the  monopulse  error  curve  from  the  adapted 
sum  and  difference  beams  in  which  the  interference  signals  have  been  suppressed.  The  resultant  dis¬ 
torted  error  curve  across  the  mainbeam  tracking  angle  region  was  used  directly  for  track  correction. 
This  approach  can  be  easily  implemented  and  is  used  here  to  estimate  the  angle  of  arrival  in  the  pres¬ 
ence  of  both  internal  noise  and  external  interference  sources  located  within  the  mainbeam  of  the  array 
antenna.  Monte  Carlo  simulation  is  performed  to  evaluate  the  approach,  and  the  angle  estimation 
errors  are  obtained.  These  errors  are  compared  with  the  Cramer-Rao  bound,  which  provides  the 
minimum  error  for  an  unbiased  estimator.  This  bound  is  derived  for  the  general  case  consisting  of 
unknown  signal  amplitude,  phase,  and  angle  of  arrival.  The  bound  has  been  previously  derived  by 
others  [3,4]  but  under  different  conditions. 

H.  PROBLEM  FORMULATION 

The  received  baseband  complex  signal  X k  at  the  &th  element  of  an  At- element  antenna  array,  as 
illustrated  in  Fig.  1,  is  given  by 

X.  =  S  eK2r/m  ~  '^Sin0s  4-  v  r  J(2*/X)(*  -  \)d  sin  ^ 

K  ~  Zj  ric  '  nk>  (1) 

i=l 

where 

S  is  the  complex  target  echo, 

<t>s  is  the  angle  of  arrival  of  the  target, 

ri  is  the  complex  signal  of  the  ith  interference  source, 

'Pi  is  the  angle  of  arrival  of  the  ith  interference  source, 

Nj  is  the  total  number  of  interference  sources, 
nk  is  the  thermal  noise  at  the  kth  element, 
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Fig.  1  —  Array  geometry 


d  is  the  element  spacing, 

X  is  the  wavelength,  and 

j  =  V-i . 

The  signal  is  assumed  to  be  processed  in  time-sampled,  digitized  inphase  I  and  quadrature  Q  chan¬ 
nels.  Each  complex  time  sample  across  the  array  comprises  a  snapshot.  Throughout  this  study  the 
antenna  element  spacing  is  X/2.  The  objective  is  to  estimate  <f>s,  the  angle  of  arrival  of  the  desired 
signal  measured  from  the  beam  pointing  direction,  in  the  presence  of  internal  noise  and  high-power, 
high-duty-cycle  external  interference  sources.  We  assume  that  the  target  signal  occurs  in  a  single 
range  cell  whereas  the  interference  signals  appear  in  many  range  cells. 

In  general  the  angle  estimation  technique  consists  of  two  measurements.  The  first  one  is  used  to 
determine  the  monopulse  error  calibration  curve.  It  uses  the  total  interference  signal  appearing  at 
each  antenna  array  element  from  which  the  covariance  matrix  can  be  computed.  The  second  meas¬ 
urement,  which  is  used  to  determine  the  target  angle,  is  described  later.  It  uses  the  total  signal 
received  in  the  presence  of  both  target  and  interference  sources. 

In  the  first  measurement,  the  covariance  matrix  R  is  estimated  from  many  range  cells  of  data  in 
the  absence  of  a  target  signal.  Here  we  first  assume  that  R  is  known  exactly,  which  is  equivalent  to 
averaging  over  an  infinite  number  of  snapshots.  In  Section  V  we  discuss  the  effect  of  using  an  R 
estimated  from  a  finite  number  of  snapshots  on  the  angle  estimation  accuracy.  Based  on  the  known 
interference  scenario,  the  total  interference  appearing  at  the  &th  array  element  is  given  by 


Nj 

Qk  -  E  n  e 


—  l)dsin^( 


+  «*• 


(2) 
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The  /^-element  of  the  covariance  matrix  is  obtained  by  applying  the  expected  value  operator  S: 


Rpq  =  £[QP  q\\ 


(3) 


where  (•)  denotes  the  complex  conjugate  and  T  is  the  transpose  operation.  The  covariance  matrix  can 
be  explicitly  expressed  in  the  form, 


Nj 

R  =  £ 

1=1 


1  —  j(2x/\)d  sin  e  2y  (2ir/\)</  sin 

j(2x/\)d  sin  ^  1  - j(2-K/\)d  sin  fa 

i 


i 


1  0 

•  0 

0  1 

0  • 

0  • 

•  1 

(4) 


where  a f  is  the  noise  power  of  the  ith  interference  source,  and  a2n  is  the  thermal  noise  power,  which 
is  assumed  to  be  the  same  for  each  element. 


The  tracking  beams  are  based  on  selecting  an  adjacent  pair  of  orthogonal,  uniform  illumination 
beams  that  are  generated  by  a  Butler  matrix  beamformer  transformation  [2],  The  resultant  sum  beam 
weight  vector  S0  and  the  difference  beam  weight  vector  D0  at  the  output  of  the  Butler  beamformer 
are  equivalent  to  cosine  illuminations. 

For  a  vector  E  of  the  received  echoes  on  the  A-element  array,  the  adapted  sum  and  difference 
beam  outputs  obtained  by  applying  conventional  beam  weighting  S0  and  D0  to  the  spatial  filtered  out¬ 
put  residue  signal  vector,  Ef  (Ef  ~  R~lE),  are  given  by 


E  =  Si  R~l  E 
and 


(5) 


A  =  Dt0  R~x  E.  (6) 

As  a  result  of  the  Butler  matrix  transformation,  A  is  90°  out  of  phase  with  E  and  the  desired  mono¬ 
pulse  ratio  is  the  imaginary  part  of  (A/E),  denoted  by  Im(A/E).  For  calibration,  Im(A/E)  is  calcu¬ 
lated  by  letting  the  vector  E  be  a  steering  vector  steered  across  the  beam  to  give  the  distorted  mono¬ 
pulse  curve  in  the  tracking  angle  region.  If  there  is  no  external  interference  source,  there  is  no  dis¬ 
tortion.  The  difference  pattern  has  a  null  at  the  beam  center  where  the  sum  pattern  peaks  (Figs.  2(a) 
and  2(b)).  The  monopulse  error  curve  is  ideally  linear  over  the  tracking  region  within  the  3  dB 
beamwidth  ($bw)>  as  shown  in  Fig.  2(c)  for  an  8-element  array. 
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Fig.  2  —  No  interference:  (a)  sum  pattern;  (b)  difference 
pattern;  (c)  monopulse  error  curve 


The  calibration  curves  for  estimating  the  angle  of  arrival  in  the  presence  of  external  interference 
are  obtained  as  follows.  Assume  that  E  is  the  desired  steering  vector, 


E  = 


1 

j{ 2x/\)  d  sin </>f 


(7) 


K 2*/X)  (N-\)d  sin <f>s 


which  is  a  function  of  angle  of  arrival  only.  Substituting  Eq.  (7)  into  Eqs.  (5)  and  (6),  one  obtains 
the  ratio  Im(A/£)  as  a  function  of  direction  of  arrival  <t>s.  For  example,  consider  a  covariance  matrix 
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being  constructed  based  on  an  interference  source  located  at  0.5  0BW  (0Bw  =  16.8815°  for  an  8- 
element  array).  The  adapted  E  and  A  patterns  are  obtained  by  using  Eqs.  (5)  and  (6)  and  are  shown 
in  Figs.  3(a)  and  3(b).  Note  that  a  notch  is  placed  at  the  corresponding  interference  source  location 
in  both  patterns  that  are  distorted  from  the  interference-free  cases  (Figs.  2(a)  and  2(b)).  Figure  3(c) 
shows  the  ratio  Im(A/E)  as  a  function  of  4>s.  The  slope  of  the  monopulse  error  curve  is  significantly 
different  from  that  shown  in  Fig.  2(c).  It  has  a  singular  point  at  the  interference  source  location.  In 
the  two-source  case  where  an  additional  interference  source  is  placed  at  —0.5  0B w>  the  adapted  E  and 
A  patterns  are  depicted  in  Figs.  4(a)  and  4(b),  respectively.  Note  that  both  patterns  have  nulls  at 
these  two  source  locations.  The  monopulse  error  curve  shown  in  Fig.  4(c),  also  different  from  Fig. 
2(c),  has  two  singular  points  where  the  interference  sources  are  located.  The  distorted  monopulse 
error  curves  computed  as  described  above  are  used  as  the  calibration  curve  in  estimating  the  angle  of 
arrival  when  the  received  signal  is  contaminated  by  the  interference  sources. 


Fig.  3  —  One  interference  source  at  0.5  beamwidth:  (a) 
sum  pattern;  (b)  difference  pattern;  (c)  monopulse  error 
curve 
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Fig.  4  —  Two  interference  sources  at  0.5  and  —0.5 
beamwidth,  respectively:  (a)  sum  pattern;  (b)  difference  pat¬ 
tern;  (c)  monopulse  error  curve 


The  second  measurement  involves  the  received  signals  consisting  of  the  returns  from  the 
interference  sources  and  the  desired  signal.  Substituting  the  vector  E,  as  given  by 

E  =  [XxX2..  .Xn]t. 


where  Xk  is  defined  in  Eq.  (1),  into  Eqs.  (5)  and  (6),  one  can  calculate  the  distorted  ratio  Im(A/E). 
The  corresponding  angle  of  arrival  <t>s  in  the  interference  environment  can  be  estimated  from  the  ratio 
Im(A/E)  based  on  the  calibration  curve  obtained  in  the  first  measurement  described  above. 

m.  THE  CRAMER-RAO  BOUND  ON  ANGLE  ESTIMATION  ERROR 

The  Cramer-Rao  (C-R)  bound  is  the  lower  bound  on  the  variance  of  any  unbiased  estimator.  It 
is  useful  to  compare  the  performance  of  an  estimation  procedure  with  this  bound.  Brennan  [3],  using 


6 


NRL  REPORT  9234 


the  C-R  bound  on  angle  estimation  error,  has  determined  the  limit  on  angle  estimation  accuracy  for  a 
target  in  the  presence  of  internal  noise  with  known  signal  amplitude  but  unknown  phase.  For  a  target 
with  an  unknown  amplitude  that  is  Rayleigh  distributed,  and  an  unknown  phase  that  is  uniformly  dis¬ 
tributed  between  0  and  2tc,  McGarty  [4]  has  derived  the  C-R  bound  in  the  presence  of  more  than  one 
interference  source  based  on  the  specific  interference  model  that  was  considered  in  obtaining  the 
covariance  matrix. 

A  general  case  consisting  of  unknown  signal  amplitude,  phase,  and  angles  of  arrival  is  con¬ 
sidered  here.  The  problem  of  angle  estimation  is  similar  to  that  described  in  Ref.  5  in  terms  of  a 
Doppler  velocity  estimation  process.  In  the  Appendix  we  show  that  for  K  independent  samples  of 
Im(A/E)  the  C-R  bound  for  the  estimation  error  of  sin  <j>s  is  given  by  (Eq.  (A10)) 


where 


(<v)2  - 


1 

SK 


cc 

CC(D  +  D)  —  AA(C  +  C ) 


n'  =  sin  <f>s. 


\yT  =  [1  ej(2*/\)d  ej(lT/\)ld  .  '  '  ej(2TT/\)(N-l)d j 

c  =  WT  R~l  W, 

A  =  WT  B  R~l  W, 

D  =  WT  B  R~l  B  W, 


0 

0 

• 

j(2Tr/\)d 

0 

• 

• 

• 

j(2ir/\)2d  0 

0 

0 

• 

0 

0 

0 

• 

•  j(2w/\)(N  -  1  )d 

(8) 


S  is  the  input  signal  power  in  each  antenna  element,  and 


K  is  the  number  of  independent  samples  of  Im  (A/E)  processed. 

This  expression  is  for  an  array  of  N  identical  antenna  elements  with  equal  spacing  d. 

A.  Single  Interference  Source: 

The  C-R  bound  on  rms  sin  <t>s  estimate  error  normalized  to  sin  0BW,  or  cy /sin  dBW,  which  is 
denoted  by  a is  calculated  as  a  function  of  angular  separation  between  the  target  and  the  interference 
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source  normalized  to  SBW.  The  scenario  consists  of  a  target  with  10  dB  signal-to-noise  ratio  (S/N0  ) 
and  an  interference  source  with  30  dB  interference-to-noise  ratio  (I/N0  ).  Figure  5(a)  shows  the  C-R 
bound  for  an  8-element  array  when  four  different  values  of  K  (K  =  1,8,  16,  and  32)  are  processed. 
It  demonstrates  that  higher  precision  in  angle  estimation  accuracy  is  achieved  as  more  samples  are 
averaged  before  estimating  the  angle  of  arrival.  The  angle  estimation  error  varies  very  little  when  the 
angular  separation  between  the  interference  source  and  the  target  becomes  slightly  greater  than  one 
beamwidth.  Generally  this  error  increases  as  the  interference  source  approaches  the  target  direction. 
The  C-R  bound  on  ,  which  is  defined  as  the  angle  estimation  error  in  beamwidths,  is  shown  in 
Fig.  5(b)  to  be  in  excellent  agreement  with  Fig.  5(a)  for  the  same  scenario.  The  vertical  axis  in  all 
the  plots  that  appear  later  in  terms  of  <rM  can  be  viewed  approximately  as  the  angle  estimation  error  in 
beamwidths. 

Figure  6  shows  the  results  for  three  different  array  sizes  when  32  independent  samples  are  pro¬ 
cessed.  The  variation  in  angle  estimation  error  with  the  angular  separation  is  similar  for  different 
array  sizes.  Also,  the  error  incurred  in  the  angle  estimation  is  less  for  larger  array  size.  Figure  7 
shows  the  C-R  bound  on  of  a  10  dB  S/N0  target  using  an  8-element  array  as  a  function  of  angular 
separation  from  an  interference  source  with  I/N0  of  10,  20,  30,  40,  and  50  dB  based  on  32  indepen¬ 
dent  samples.  It  shows  that  the  angle  estimation  error  becomes  more  severe  as  the  power  of  the 
interference  source  increases.  For  an  angular  separation  of  more  than  approximately  one  beamwidth, 
the  rms  error  levels  off  and  the  power  of  the  interference  source  does  not  significantly  affect  the  esti¬ 
mation  error. 

The  scenario  in  Fig.  8  is  identical  to  that  in  Fig.  7  except  the  target  has  a  S/N0  of  0  dB  and  the 
interference  source  has  four  different  values  of  I/N0  (I/N0  =  0,  10,  20,  and  30  dB).  Again  with  an 
angular  separation  of  slightly  greater  than  one  beamwidth,  the  estimation  error  is  the  same  for  all 
I/N0  power  levels.  A  comparison  of  Fig.  8  with  Fig.  7  shows  that  better  angle  estimation  is  achieved 
for  higher  S/N0  with  the  same  I/S  ratio  in  both  cases.  Figure  9  is  the  same  as  Fig.  8  except  N  =  16. 
It  is  demonstrated  by  comparing  these  two  figures  that  a  smaller  error  is  obtained  for  a  larger  array 
size  with  the  identical  target  and  interference  source  scenarios. 

B.  Two  Interference  Sources: 

In  these  calculations,  it  is  assumed  that  the  target  is  at  the  center  of  the  antenna  mainbeam.  The 
scenario  considered  here  is  that  the  two  interference  sources  of  equal  power  are  located  symmetrically 
with  respect  to  the  target.  Figure  10  shows  the  C-R  bound  of  plotted  as  a  function  of  angular 
separation  between  the  target  and  one  of  the  two  symmetrically  located  interference  sources,  normal¬ 
ized  to  6bw,  for  four  different  values  of  A"  (A"  =  1,8,  16,  and  32).  The  parameters  used  in  Fig.  10 
are:  N  =  8,  S/N0  =  10  dB,  and  I/N0  =  30  dB  for  both  interference  sources.  The  error  is  reduced 
as  the  interference  sources  simultaneously  move  away  from  the  target  (from  0.1  to  1.0  dBW).  The 
error  approaches  a  constant  value  as  the  interference  sources  are  more  than  one  beamwidth  from  the 
target. 

IV.  ANGLE  ESTIMATION  ERROR  BASED  ON  SIMULATION  RESULTS 

The  method  described  in  Section  II  is  used  to  estimate  the  angle  of  arrival  in  the  presence  of 
mainbeam  interference.  The  normalized  angular  error  is  then  compared  with  the  bound.  The 
received  signal  Xk  as  given  in  Eq.  (1)  is  simulated  for  each  sample  of  a  run.  Random  complex  Gaus¬ 
sian  variables  are  generated  to  represent  interference  sources  and  receiver  noise.  The  E  and  A  beam 
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ANGULAR  SEPARATION  IN  BEAMWIDTHS 

(a) 


(b) 

Fig.  5  —  (a)  The  Cramer-Rao  bound  of  ah  ,  (b)  The  Cramer-Rao  bound  of 
oc, ,  as  a  function  of  angular  separation  for  different  values  of  K:  N  =  8; 
S/N0  =  10  dB;  I/N0  =  30  dB. 
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Fig.  6  —  The  Cramer-Rao  bound  of  a M  as  a  function  of  angular  separa- 
ation  for  different  values  of  N:  K  =  32;  S/N0  =  10  dB;  I/N0  =  30  dB. 


Fig.  7  —  The  Cramer-Rao  bound  of  o ^  as  a  function  of 
angular  separation  for  different  values  of  I/N0:  TV  =  8;  K 
=  32;  S/N0  =  10  dB. 
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ANGULAR  SEPARATION  IN  BEAMWIDTHS 

Fig.  8  —  The  Cramer- Rao  bound  of  a rM  as  a  function  of  angular 
separation  for  different  values  of  I/N0:  N  -  8  ;  K  =  32;  S/N0  = 
OdB. 


Fig.  9  —  The  Cramer-Rao  bound  of  <rM  as  a  function  of 
angular  separation  for  different  values  of  I/N0:  N  =  16;  K 
=  32;  S/N0  =  0  dB. 
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Fig.  10  —  The  Cramer- Rao  bound  of  o M  as  a  function  of  angular 
separation  for  different  values  of  K  in  the  presence  of  two  interference 
sources  with  target  at  the  beam  center:  N  =  8;  S/N0  =  10  dB;  Ij/No  — 
30  dB;  I2/N0  =  30  dB. 


returns  are  obtained  by  using  Eqs.  (5)  and  (6)  for  each  run,  and  the  ratio  Im(A/E)  is  calculated. 
Finally,  ^-independent  samples  of  the  ratio  are  averaged  before  the  target  angle  is  estimated  by  using 
the  calibration  curves  obtained  from  the  first  procedure  of  the  estimation  technique  described  in  Sec¬ 
tion  II.  The  results  presented  here  are  based  on  repeating  the  above  process  of  angle  estimation  100 
times. 

In  the  single-source  case  the  following  scenario  is  assumed.  An  antenna  array  with  8  elements 
is  used  to  receive  the  signal  consisting  of  the  returned  echo  from  a  target  with  S/N0  =  10  dB  and  an 
interference  source  with  I/N0  =  30  dB.  The  interference  source  location  is  chosen  first  and  the  tar¬ 
get  direction  is  varied  within  the  mainbeam.  Figures  11(a),  11(b),  11(c),  and  11(d),  respectively, 
show  the  results  for  the  selected  interference  source  locations  (0.2,  0.5,  0.7,  and  0.9  dBW,  for  exam¬ 
ple)  as  the  target  direction  is  varied  from  the  interference  source,  first  toward  the  beam  center  and 
then  away  from  it.  The  solid  lines  are  the  theoretical  C-R  bounds  on  a and  the  individual  data 
points  are  from  the  simulation.  Four  different  values  of  K  (K  =  1,  8,  16,  and  32)  are  used  to 
demonstrate  the  effects  on  the  angle  estimation  accuracy.  All  the  simulated  points  are  near  the  C-R 
bounds  for  all  values  of  K  tested.  Figures  12(a),  12(b),  12(c),  and  12(d),  for  the  four  different  values 
of  K,  respectively,  show  the  combined  results  with  the  interference  source  locations  varying  from  0.1 
to  0.9  0BW  with  a  step  of  0.1  dBW.  The  simulated  points  of  the  same  symbol  correspond  to  the 
scenario  with  a  fixed  interference  source  location  and  various  target  locations.  Here  the  target  is 
located  in  the  mainbeam  but  not  necessarily  at  the  mainbeam  center.  Figure  13  shows  the  results  for 
the  case  when  the  target  is  assumed  at  the  mainbeam  center  and  the  interference  source  location  is 
varied  from  0.2  to  0.9  dBW  with  an  increment  of  0.1  6BW  for  four  different  values  of  K  (K  =  1,8, 
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(a) 


Fig.  11  —  The  Cramer-Rao  bound  ( - )  and  the  simulated  data  points: 

K  =  1  (xxx),  K  =  8  (OOO).  K  =  16  (V  V  V),  and  K  =  32 
(G3EIKI),  of  as  a  function  of  angular  separation:  N  —  8;  S/N0  =  10 
dB;  I/N0  =  30  dB.  Interference  source  is  at:  (a)  0.2,  (b)  0.5,  (c)  0.7, 
and  (d)  0.9  beam  width. 
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K=  1 


K  =  8 
K  =  16 
K  =  32 


Fig.  11  (Continued)  —  The  Cramer- Rao  bound  ( - )  and  the  simulated 

data  points:  K  =  1  (xxx),  K  =  8  «>  O  O),  K  =  16  (V  V  V),  and  K 
=  32  (H  (E^),  of  0^  as  a  function  of  angular  separation:  N  =  8;  S/N0 
=  10  dB;  I/N0  =  30  dB.  Interference  source  is  at:  (a)  0.2,  (b)  0.5,  (c) 
0.7,  and  (d)  0.9  beamwidth. 
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ANGULAR  SEPARATION  IN  BEAMWIDTHS 


(a) 


Fig.  12  —  The  Cramer-Rao  bound  ( - )  and  the  simulated  data 

points  for:  (a)  K  =  1,  (b)  AT  =  8,  (c)  K  =  16,  and  (d)  K  =  32, 
of  as  a  function  of  angular  separation:  N  =  8  ;  S/N0  =  10  dB; 
I/No  =  30  dB.  Interference  source  location  is  varied  from  0.1  to 
0.9  beamwidth. 
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(d) 

Fig.  12  (Continued)  —  The  Cramer-Rao  bound  ( - )  and  the 

simulated  data  points  for:  (a)  K  —  1,  (b)  K  =  8,  (c)  K  —  16,  and 
(d)  K  =  32,  of  0^  as  a  function  of  angular  separation:  N  =  8  ; 
S/N0  =  10  dB;  I/N0  =  30  dB.  Interference  source  location  is 
varied  from  0.1  to  0.9  beamwidth. 
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ANGULAR  SEPARATION  IN  BEAMWIDTHS 

Fig.  13  —  The  Cramer-Rao  bound  ( - )  and  the  simulated  data  points: 

K  =  1  (xxx),  K  =  8  (OOO).  K  =  16  ( V  V  V ),  and  K  =  32 
(53EE3IS3),  of  as  a  function  of  angular  separation  with  target  at  the  beam 
center:  N  =  8;  S/N0  =  10  dB;  I/N0  =  30  dB. 


16,  and  32).  The  C-R  bound  of  the  angle  estimation  error  is  also  plotted  in  Fig.  13.  The  angular 
error  obtained  from  the  estimation  procedure  described  here  for  the  single-source  case  is  shown  to  be 
quite  close  to  the  C-R  bound  derived  here.  The  error  is  less  than  0. 1  dBW  for  K  >  8. 

Simulation  is  also  performed  for  the  two-source  case.  The  scenario  considered  here  includes  a 
target  at  the  mainbeam  center  and  two  interference  sources  symmetrically  located  on  both  sides  of  the 
target  direction  within  the  mainbeam.  The  parameters  used  here  are  the  same  as  in  the  single-source 
case,  except  there  are  two  interference  sources  with  equal  I/N0  of  30  dB.  Figure  14  shows  the  simu¬ 
lation  results  and  the  C-R  bounds.  For  an  angular  separation  between  the  target  and  one  of  the  sym¬ 
metrically  located  interference  sources  greater  than  0.5  0B\ y,  the  simulated  points  obtained  by  the  esti¬ 
mation  method  described  here  are  within  0.1  6BW  for  K  >  8  and  are  close  to  the  C-R  bound.  As  the 
separation  between  the  target  and  the  interference  sources  becomes  equal  to  or  less  than  0.5  dBW,  the 
accuracy  of  angle  estimation  starts  to  degrade. 

V.  EFFECTS  OF  ESTIMATING  THE  COVARIANCE  MATRIX 

So  far  we  have  assumed  that  the  covariance  matrix  is  obtained  by  averaging  over  an  infinite 
number  of  range  cells.  In  practice,  the  covariance  matrix  must  be  estimated  from  a  finite  number  of 
snapshots.  The  effect  of  the  accuracy  of  the  estimated  covariance  matrix  on  the  accuracy  of  the  angle 
estimation  is  briefly  examined  here.  & 
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Fig.  14  —  The  Cramer-Rao  bound  ( - )  and  the  simulated  data 

points:  K  =  1  (xxx),  K  =  8  (<0>OO)»  K  =  16  (V  V  V), 
and  K  =  32  (H  81  Kl),  of  as  a  function  of  angular  separation 
in  the  presence  of  two  interference  sources  with  target  at  the 
beam  center:  N  =  8;  S/N0  =  10  dB;  I,/N0  =  30  dB;  I2/N0  = 
30  dB. 


First,  the  covariance  matrix  R  is  estimated  by  using  32  snapshots.  Procedures  for  estimating  the 
angle  of  arrival  of  the  desired  signal  previously  described  are  repeated  by  using  the  estimated  rather 
than  the  ideal  value  of  R.  Figure  15(a)  shows  the  results  obtained  for  an  interference  source  located 
at  0.5  6bw  with  the  same  scenario  as  that  used  in  Fig.  11(b).  Solid  lines  are  the  C-R  bounds,  and  the 
individual  data  points  are  from  the  simulation.  The  simulated  angle  errors  are  far  above  the  C-R 
bounds. 

Next,  256  snapshots  are  used  for  estimating  R.  Figure  15(b)  shows  the  results.  Note  that  the 
errors  in  angle  estimation  are  much  lower  than  those  shown  in  Fig.  15(a).  As  the  number  of 
snapshots  increase,  the  angle  estimation  accuracy  improves  and  approaches  what  is  shown  in  Fig. 
11(b)  for  an  idealized  covariance  matrix.  Also,  more  snapshots  are  required  to  approach  the  bound  as 
the  angular  separation  decreases. 

VI.  SUMMARY 

An  investigation  has  been  performed  on  estimating  the  angle  of  arrival  of  a  desired  signal  in  the 
presence  of  mainbeam  interference.  Adaptive  antenna  arrays  are  incorporated  to  form  adapted  sum 
and  difference  beams  in  which  the  interference  signals  are  suppressed.  Monopulse  error  curves  are 
then  obtained,  providing  the  necessary  distortion  correction  curves  across  the  entire  mainbeam  track¬ 
ing  angle  region.  New  C-R  bounds  on  the  angle  estimation  error  are  derived  with  generalized 
assumptions  on  the  signal  amplitude  and  phase.  The  bounds  previously  derived  by  others  are  valid 
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Fig  15  The  Cramer-Rao  bound  ( - )  and  the  simulated  data  points: 

K  =  1  (xxx),  K  =  8  (OOO).  ^  =  16  (V  V  V)  and  K  =  32 
(EIKIKI),  of  as  a  function  of  angular  separation  for  an  interference 
source  located  at  0.5  beamwidth:  N  =  8;  S/N0  =  10  dB;  I/N0  =  30  dB. 
R  is  estimated  from  (a)  32;  (b)  256  samples. 
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under  different  conditions.  With  these  generalized  assumptions  on  the  signal  characteristics,  a  Monte 
Carlo  simulation  is  performed  based  on  the  estimation  procedure  described  herein  to  determine  the 
angle  estimation  error.  These  errors  are  compared  with  the  C-R  bounds.  Good  performance  is 
shown  for  sufficient  S/N0  and  angular  separation  between  the  target  and  the  interference  sources. 

For  the  case  of  an  8-element  array  with  a  single  interference  source,  the  angle  estimation  error 
is  within  0.1  dBW  when  the  angular  separation  is  equal  or  greater  than  0.2  dBW  and  at  least  eight 
independent  samples  of  the  monopulse  measurement  Im(A/E)  are  processed.  To  achieve  the  same 
accuracy  in  angle  estimation  for  the  two-source  case,  the  angular  separation  must  be  greater  than  0.5 

®BW- 


The  issue  of  the  effect  of  estimating  the  covariance  matrix  from  a  finite  number  of  snapshots  on 
the  angle  estimation  accuracy  is  also  briefly  addressed.  It  is  shown  that  the  accuracy  in  angle  estima¬ 
tion  is  greatly  affected  by  the  estimated  covariance  matrix.  In  general,  more  snapshots  are  needed  to 
approach  the  bounds  as  the  angular  separation  between  the  target  and  the  interference  source 
decreases. 
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Appendix 


Let  X,-  be  an  A-dimensional  complex  signal  vector  representing  the  z'th  sample  of  the  signal 
received  by  the  iV-element  linear  antenna  array.  The  received  signal  vector  X, ,  composed  of  a  signal 
component  and  a  noise  component  n, ,  is  given  by 

%  =  st  +  rti  =  b  W  +  rij,  (Al) 

where  /z,  is  a  complex  zero-mean  Gaussian  process  with  covariance  matrix  R,  b  is  the  unknown  com¬ 
plex  signal  strength,  and  W  is  the  steering  vector.  For  an  iV-element  linear  antenna  array  of  identical 
elements  and  uniform  spacing  d,  the  steering  vector  is  written  in  terms  of  A,  the  phase  shift  from  ele¬ 
ment  to  element,  owing  to  the  angle  of  arrival  <t>s : 

WT  =  [1  ejA  e2JA  . . .  e(N~1)jA],  (A2) 

where  T  denotes  the  transpose,  j  is  the  square  root  of  -1,  A  =  (2ir/\)d  sin  <f>s,  and  X  is  the 
wavelength.  The  vector  X,  has  a  conditional  probability  density 

P  (X,  I  ©)  =  lR  exp  [  -  (X,  -  bW)  R  -1  (X,-  -  bW)},  (A3) 


where  0  represents  the  unknown  parameters  =  sin  b,  and  b;  (•)  denotes  the  complex  con¬ 
jugate;  and  |  R  |  is  the  determinant  of  the  noise  covariance  matrix.  The  joint  probability  density  of  K 
independent  samples  can  be  written  as 


p  (Xj  X2,  ...  ,  X*|0)  =  — 1  exp  [  -  £  (X,  -  bW)T  R-'  (X;  -  bW)}.  (A4) 

*  I"  I  i=i 

The  Cramer-Rao  bound  on  the  variance  of  /z'  is  obtained  from  the  Fisher’s  information  matrix  J 
whose  elements  are 


Ju  =  E 


3  In  PiXu 


O) 


36 


d  In  j7(Xlt  ...  Xg|  0) 
ddj 


(A5) 


where  6i  —  p',  02  =  b,  03  =  b\  E  is  the  expected  value  operator;  and  i  and  j  are  indices.  By  using 
the  density  function  (Eq.  (A4))  and  the  definition  (Eq.  (A5)),  then 
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Note  that 


J  ll 


=  bb 


'dW  i  dW_  dW  ! 

dp'  dp'  dp'  dp' 


K 


J2 2  =  K  WT  R~x  W  =  K  C 
J33  =  K  WT  R~x  W  =  K  C 

Jn  =b  K  4tT"  *-1w  =  b  KA 

Ofl 

—  T 

jn  =bK^y-rR~xW  =  bKA 

J2 1  =  b  K  WT  R~x  4^7  =  & 

•  dp 

J3l  =  b  KWT  R~x  ^-  =  b  K  A 

OfX 

J  23  =  ^32  =  0. 


air 

dp' 


=  BW  and 


dp' 


=  !Fr  B,  where  B  is  defined  as  the  iV  X  A/  diagonal  matrix: 


0  0.  0 

0  j(2ir/\)d 

2j(2ir/\)d 

0 

0  •  (N  -  l)j(2ir/\)d 


(A6) 
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The  matrix  J  can  then  be  inverted  in  a  closed  form 


r-l  - 


CC 


-b  AC 


-b  A  C 


-b  AC  CJn  -bbAA  b2  A2 


-b  AC 


b2  A2 


C  Jn  -bbAA 


(A7) 


where  |  J  \ ,  the  determinant  of  the  matrix  J,  is  expressed  by 

I J  I  =CCJn-bbAAC-bbAAC.  (A8) 


The  bound  on  the  variance  of  the  estimate  of  /z '  required  by  the  Cramer-Rao  bound  is 


a 


2 


> 


c  c 
\J  I  ’ 


(A9) 


which  is  the  upper  left-hand  comer  of  the  matrix  given  in  Eq.  (A7).  Equation  (A9)  can  be  explicitly 
written  as 

g2  > _ _ (WlrAV)(W^R  ‘wq _ _ _ 

h  5/sr[(tT7'R“1H0(tT7R“1W0(^^ ,SR_I  BtT  +  tP7'BR'l5H0  -  +  Wlc~AV)]  ' 

(A10) 


where  S  =  bb  is  the  signal  power. 
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